Nonextensive Thermostatistics and the 
//-Theorem Revisited 



Fernando M. Ramos Reinaldo R. Rosa* 
Luis A. W. Bambace* 

^Instituto Nacional de Pesquisas Espaciais, INPE 
12201-970 Sao Jose dos Campos - SP, Brazil 



Abstract 

In this paper we present a new derivation of the H-theoiem and the correspond- 
ing colhsional equiUbrium velocity distributions, within the framework of Tsahis' 
nonextensive thermostatistics. Unhke previous works, in our derivation we do not 
assume any modification on the functional form of Boltzmann's original " molecular 
chaos hypothesis". Rather, we explicitly introduce into the collision scenario, the 
existence of statistical dependence between the molecules before the collision has 
taken place, through a conditional distribution f{v2\vi). In this approach, different 
equilibrium scenarios emerge depending on the value of the nonextensive entropic 
parameter. 



Boltzmann's kinetic theory, certainly one of the most important developments 
in the history of statistical mechanics, has as one of its key ingredients the 
"molecular chaos hypothesis", the celebrated (and somewhat controversial) 
Stosszahlansatz. The irreversible nature of Boltzmann's equation is directly 
linked to this hypothesis, since molecules in a dilute gas are assumed to be 
correlated only after the collision has taken place [1] . 

Recently, Lima et al. [2] (see also [3]) advanced an ad-hoc modification of the 
Stosszahlansatz, in order to derive a proof of the if-theorem, and to obtain the 
corresponding colhsional equilibrium velocity distribution, within the frame- 
work of Tsallis' nonextensive thermostatistics [4] . In this paper, we show that 
there is no need of modifying the functional form of the molecular chaos hy- 
pothesis for accomplishing this task. All we need, besides adopting Tsallis' 
prescription for the functional form of the local entropy, is to explicitly as- 
sume that molecules are also correlated before the collision. In other words, 
we assume that the joint distribution /(I, 2) associated with two colliding 
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molecules, factorizes as /(1, 2) = /i^2 = /(1)/(2|1) = /1/21, where /21 con- 
tains all the information regarding the correlations between 1 and 2, before 
the impact. 

Let us start considering a spatially homogeneous gas of N hard-sphere particles 
of mass m and diameter s, in the absence of an external force. The state of 
a dilute gas is kinetically characterized by the velocity distribution function 
f{v,t). The quantity f{vi,V2,t)(fvicPv2 gives, at each time t, the number of 
particles in the volume element (Pvi around the velocity vi and in the volume 
element (fv2 around the velocity V2- The rate of change of f{vi,V2,t)(Pvid^V2 
verifies the generalized Boltzmann equation 

% = Wl.2) , (1) 



where Cq denotes the generalized collisional term. We also make the usual 
assumptions: (i) Only binary collisions occur in the gas; (ii) Cg(/) is a local 
function of the slow varying distribution function; (iii) Cq{f) is consistent with 
the energy, momentum, and particle number conservation laws. Now, following 
standard lines we define 
2 

Cq{fi,2) = j\V-e\Rq du;d''v'id'v'2 , (2) 



where the primes indicate post collision quantities, 1/ = 1/2 — tTi is the relative 
velocity before collision, e denotes an arbitrary imit vector, and doj is an 
elementary solid angle such that s'^du) is the area of the collision cylinder 
[5,6]. The quantity Rq{f,f') is a difference of two correlation functions (just 
before and after collision). At this point we depart from [2], who originally 
proposed that /') = exp,(/(^-^ In, fi) + fr' K /^) -exp,(/r^ In, f,) + 
/|~^ In, /2), where exp^ and In, are ^-exponential and ^-logarithmic functions 
(see below), and assume that /?,(/, /') simply satisfies the stanc/aro? functional 
form of Boltzmann 's molecular chaos hypothesis 

W,/')=/(/21-/l/21 , (3) 



where again primes refer to the distribution function after collision. Note that, 
although the functional form is the same, there is a subtle but profound differ- 
ence between ours and Boltzmann's original approach, in the sense that now 
velocities of colliding molecules are not statistically independent before the 
impact. 

For the joint entropy we adopt Tsallis' expression, 

H'/ ^-kj fl2 In, /i,2 d'v,d'v2 , (4) 



2 



where liio = -^^ ^ ^ . Note that Hl''^ reduces to the standard Boltzmann exten- 

1 l—q q 

sive measure if^'^ — + H'^ , ior q — 1. Taking the partial time derivative of 
the above expression we obtain 



dt 



-k J{qf!;,'ln,h,2 + l)^d'v,d'v2 . (5) 



Inserting the generalized Boltzmann equation (1) into (5), and using (2), ex- 
pression (5) can be rewritten as a balance equation 

-^ = Gq{t) , (6) 



where the source term Gg reads 

Gg = -— j\V-e\{l + qfl^' In, h^^)Rg du^v.^v^d^'^dh'^ . (7) 



In order to rewrite Gq in a more symmetrical form, some elementary operations 
must be done in the above expression. Following standard lines [5], we first no- 
tice that interchanging vi and V2 does not affect the value of the integral. This 
happens because the magnitude of the relative velocity vector and the scatter- 
ing cross section are invariants. Similarly, we may use time-reversal symmetry. 
Note that this step requires the change of sign of (inverse collision). Imple- 
menting these operations and symmetrizing the resulting expression, one may 
show that the source term can be written as 

Gq^-^l \V-^{qf!:2'Kfi,2-qfS2'Kf[,2) 

Rqdu;d%d%d^v'id^v'2 • (8) 

Noting that /i^2 = /1/21, and applying the transformation /^^^ In, / = In,* /, 
where q* — 2 — q, we finally find 

Gg^^JlV- e\{\nMf^,) - WifMMfL - /1/21) 

dujd^vid^V2d^v\d^v'2 ■ (9) 



For positive values of q, the g-log function is concave. Thus, the integrand in 
equation (9) is never negative because (In,, /(/^i-ln,. /1/21) and (/(/2i-/i/2i) 
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always have the same sign. Therefore, we obtain the iJ^-theorem 



dH. 



1,2 



dt 



>0 , (10) 



for g* > or, equivalently, < g < 2. This inequality describes the route 
to equilibrium within the framework of Tsallis' nonextensive formalism. We 
remark that, since /i and /2 arc not statistically independent, we cannot say 
anything about (or H^) individually, whose value may even decrease during 
the relaxation process. 

Now, to obtain the equilibrium g-distribution all we need is to set Gq — 0. 
Since the integrand appearing in (9) cannot be negative, the source term will 
vanish, for any pair of molecules obeying the energy, momentum, and particle 
number conservation laws, if and only if 

/;/^i = /i/2i , (11) 

or equivalently 

lUq f[ + lUg /21 + (1 - g) In^ f[ \nq /21 = \nq fi + lUq /21 + 

(l-g)ln,Aln,/2i . (12) 



Clearly, In^ fi and In^ /21 + (1 — ?) In^ fi In^ /21 arc collision invariants. Thus, 
the marginal equilibrium distribution must be of the form 

In^ fi^ ao + ai-vi + a2Vi ■ vi , (13) 



where Oq and 02 are constants and ai is an arbitrary constant vector. By 
introducing the bary centric velocity u, we may rewrite (13) as 

In^ fi = a- ^\vi - , (14) 



with a different set of constants. Now, defining Tq — eq{a) and = (jqir^fi^)^ ; 
we obtain a g-generalized Maxwell's distribution 

foivi) = r, eq{-P\v, - u\') = r, [1 - (1 - q)P\v, - ul']'/'"! , (15) 



where eq{x) = [1 + (1 — q)xY^^ and Tq, (3 and u may be functions of the 
temperature. This result is identical to that obtained in Ref. [2], using a very 
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different functional form for the molecular chaos hypothesis. Naturally, Boltz- 
mann's classical scenario, and its corresponding results, are recovered in the 
limit of g — > 1. 

Now, assuming that, by collision invariance, we have in equilibrium In^ /21 + 
(1 — g) lug fi lug /21 ~ —l3{\v2 — the g-generalized conditional and joint 
distributions can be easily obtained 



and 



+ \V2 - U\^)] 



(17) 



with Qq and Ag being the appropriate normalization factors. 

Note that, although Vi and V2 are, by definition, statistically dependent^ is 
easy to verify that the covariance {vi,V2) equals zero, regardless the value of 
q. A natural question arises from this somewhat paradoxycal result: what is 
therefore the effect induced by nonextensivity on the classical picture of a 
dilute gas? 

To answer this question we plotted in Fig. 1 the function Xq{'^ij''^2), where 
/i,2 = X9/1/2, for three values of the entropic parameter. Xq ^cts as a g- 
generalized, velocity-dependent Enskog factor, reducing or enhancing the col- 
lision frequency at certain velocities. For g > 1, high- velocity collisions are 
favored, the rms velocity increases, the mean collision time decreases, and 
the velocity distribution ceases to be a Gaussian, displaying power-law tails. 
Strong, intermittent fluctuations on the velocity and the particle density fields 
are expected to emerge. This behavior is reminiscent to that observed in gran- 
ular gases [7]. For q < 1, we have the opposite scenario: high- velocity collisions 
arc filtered, the rms velocity decreases, the mean collision time increases, and 
the velocity distribution displays a cut-off, meaning that the gas occupies 
only a subset of its phase space. A condensate of cold, motionless particles 
may eventually appear. Naturally, for g = 1, we have Xq — ^, regardless the 
velocities, what recovers Boltzmann's original molecular chaos hypothesis. 

Summarizing, in this paper we presented a new derivation of the if-theorem 
and the corresponding coUisional equilibrium velocity distributions, within the 
framework of Tsallis' nonextensive thermostatistics. Unlike previous works, in 
our derivation we do not assume any modification on the functional form 
of Boltzmann's original "molecular chaos hypothesis". Rather, we explicitly 
introduced into the collision scenario, the existence of statistical dependence 
between the molecules before the collision, through a conditional distribution 
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Fig. 1. Variation of Xq for = 1-2 (top), q = 1 (intermediate) and q = 0.8 (bottom). 

fiv2\vi). In this approach, different equihbrium scenarios emerge depending 
on the value of the nonextensive parameter. 
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